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ABSTRACT 

We present a model for generation mechanisms of the optical pulsations recently de- 
tected from the Geminga pulsar. We argue that this is just a synchrotron radiation 
emitted along open magnetic field lines at altitudes of a few light cylinder radii (which 
requires that Geminga is an almost aligned rotator), where charged particles acquire 
non-zero pitch-angles as a result of the cyclotron absorption of radio waves in the 
magnetized pair plasma (Gil, Khechinashvili & Melikidze 1998, hereafter Paper I). 
This explains self-consistently both the lack of apparent radio emission, at least at 
frequencies higher than about 100 MHz, and the optical pulsations from the Geminga 
pulsar. From our model it follows that the synchrotron radiation is a maximum in 
the infrared band, which suggests that Geminga should also be a source of a pulsed 
infrared emission. 
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1 INTRODUCTION 

, An optical counterpart of the Geminga source was found by 
■ Bignami et al. (1987), based upon color considerations. Re- 
cent spectral studies (Martin et al. 1998) show a continuous 
power-law (^^-o s±o.s^ ^^.^^ 37qq qqqq ^ ,^^^^1 a broad ab- 

' sorption feature over 6300 -r- 6500 A band. Such a spectrum 
, indicates existence of the dominant nonthermal synchrotron 
' component in the optical emission of the Geminga pulsar. 

Recently Shearer et al. (1998), based on their deep inte- 
' grated images in _B-band, claimed that the optical counter- 
part of Geminga pulses with a period P — 237 ms, reported 
also at X (Halpern & Holt 1992) and 7-ray (Bertch et al. 
1992) bands. They derived the magnitude of the pulsed emis- 
sion rrig = 26.0 ± 0.4 mag, and found that the signal shows 
two peaks with a phase separation of « 0.5. It appeared that 
there is a phase agreement of the optical data with 7 and 
hard X-ray curves. Moreover, the form of the optical light 
curve resembles that of 7 and hard X-rays rather than the 
soft X-rays signature. According to Shearer et al. (1998), 
this is a further evidence of the predominantly magneto- 
spheric origin of Geminga's pulsed optical emission (Halpern 
et al. 1996) over a thermal one (Bignami et al. 1996). 

We suggest that Geminga's pulsed optical emission is 
a synchrotron radiation of the plasma particles gyrating 
about relatively weak magnetic field of the distant mag- 
netosphere of this pulsar. As we demonstrate below, the 
plasma particles acquire perpendicular momenta, necessary 
for the synchrotron emission, due to cyclotron absorption of 
radio waves on these particles. There is an evidence that the 
Geminga pulsar is either radio quiet or, more likely, visible 



only at low radio frequencies below about 100 MHz (e.g. 
Malofeev & Malov 1997; Vats et al. 1999; see also Paper I). 



2 ESTIMATION OF PITCH-ANGLES 

Mikhailovskii et al. (1982) were first to suggest that the en- 
ergy lost by radio waves due to cyclotron damping should 
lead to 'heating' of the ambient plasma, i.e. increase of 
its perpendicular temperature. Lyubarskii & Petrova (1998) 
treated thoroughly spontaneous re-emission of the absorbed 
energy, and concluded that in short-period pulsars a signifi- 
cant fraction of this energy is re-emitted in the far infrared 
band. They found that an electron (or positron), moving 
initially with the relativistic velocity strictly along the local 
magnetic field, obtains a pitch-angle^ 

Tp ^ e-n^^^ rad (1) 

in the course of cyclotron absorption of radio emission. Here 
9 is an angle between the wave vector fcr and the local mag- 
netic field B, and 

(2) 

is the ratio of the radio luminosity to the particle luminosity, 
that is, the ratio of the power of radio waves (damped due 
to cyclotron resonance) to the total kinetic energy of the 

* The pitch-angle is defined as -0 = arctan ^px/P||)i where py 
and are the components of the particle momentum along and 
across the local magnetic field, respectively. 
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Figure 1. Schematic representation of the region where the bulk 
of Geminga's radio emission is damped. This cut-cone-shaped vol- 
ume is constrained by the two dashed circles transverse to /x. The 
beaming of the optical radiation, re-emitted from this region, is 
determined by the angle (p. See definitions in the text. 



resonant particles produced per second above a fraction of 
the polar cap with the area Spi (see equation Note that 
equation (^) is vahd provided that ^ 1. Estimation of 77 
is presented below. 

The radio luminosity of the Geminga pulsar, by what we 
mean the integrated radio power of Geminga as it would be 
in the absence of cyclotron damping on the magnetospheric 
plasma (Paper I), can be estimated as 



~ F {100 MUz) Au d'^ [erg s"^], 



(3) 



where F (100 MHz) ^ IJy = 10"^^ erg cm'^ s-^Hz"^ (Vats 
et al. 1999), Au ~ 1 GHz is a frequency range where the 
bulk of pulsar radio emission is usually radiated, and d « 
160 pc « 5 X 10^" cm is a distance to the Geminga pulsar. 
Inserting all these values into equation (^ we find that Lr ~ 
3 X 10^^ erg s"^ 

Let us now evaluate the denominator in equation 
According to Paper I (see also Malov 1998), the Geminga 
pulsar is an almost aligned rotator, and for the inclina- 
tion angle a = 5° and the impact angle /3 — 20°, an ob- 
server's line of sight in fact grazes the emission cone radi- 
ated from the bundle of the last open field lines, with the 
opening angle at about 100 MHz being pi « 30° (Khechi- 
nashvili, Melikidze & Gil 1999). The latter defines a group of 
the radial lines constituting an angle pi with the magnetic 
axis fj, (see Fig. ^). The simulation of radio waves damp- 
ing along these lines (Paper I; Khechinashvili et al. 1999) 
shows that 1/ ~ 100 MHz radio waves are damped at the 
distance TZ2 ~ 2.4, whereas 1/ ^ 1 GHz waves are damped 
at TZi « 1.1. Here TZ = R/R^^ is a distance from the star's 
centre measured in the units of the light cylinder radius 
(J?Lc ~ Pc/2n). One can speculate that the major part 
of Geminga's radio emission is effectively damped within a 
cut-cone-shaped volume represented in Fig. ^. Then, let us 
consider the dipolar field line which crosses the radial line, 
constituting the angle pi — 30° with fi, at the distance 
TZi ~ 1.1. The footprint of this field line on the stellar sur- 
face has a polar angle i9pi ~ {Ro/ Ri)^'^ pi ~ 0.014 rad. The 
footprints of all inner field lines, on which 0.1 -^ 1 GHz radio 



frequencies are damped at high altitudes (see above), mark 
a circle with the angular radius i?pi within the pulsar polar 
cap through which the resonant leave the neutron star. Let 
us notice that i9pi ^ 0.5 i9p, where ~ (Rq/Rj^^)^^^ is an 
angular radius of the pulsar polar cap, and Ro ~ 10^ cm 
is a stellar radius. Hence, only a fraction of the polar cap 
provide resonant particles. The area of this region is 



; 6 X 10" cm^ 



(4) 



(5) 



Si ^ ttRo^Ii 
which yields 

Lpi = K,n^jS'yj,mc^ « 8 x 10^"^ k [erg s""'] 

Here the surface value of the corotational charge num- 
ber density (Goldereich & Julian 1969) n^, « 2.2 x 

10^* [P/Py^'^ ^ 5.4 X 10^^ cm"3 was substituted, which 
was evaluated for the dynamical parameters of the Geminga 
pulsar (P = 0.237 s and P = 1.1 x 10"^*); k is the Stur- 
rock multiplication factor (Sturrock 1971). Evaluating equa- 
tion (P) we assumed that the mean Lorentz factor of the 
ambient plasma particles 7^ ^ 100. From equation (^) it 
follows that 
^ Lpi 



"pi 



7 mc^ 



10^^ K 



(6) 



particles are generated per second over the area of Si (equa- 
tion ||). 

Substituting the derived values Lr and Lpi in equa- 
tion ^ we obtain that rj « 0.4k~^. Then, taking into ac- 
count the average value of the angle 9 ~ 0.6 rad (Khechi- 
nashvili et al. 1999) in equation (jl|), we find that the charac- 
teristic pitch-angle of plasma particles is ip ki 0.4k~'''^ rad. 
Therefore, a 'kick' that each electron gets due to the cy- 
clotron absorption of radio waves, depends on the number 
of particles (determined by the Sturrock multiplication fac- 
tor k) among which the energy of these waves is distributed. 
For example, inside the spark-associated columns of plasma 
(e.g. Ruderman & Sutherland 1975; Gil & Sendyk 1999), 
where k ~ 10*, we have ip ~ 10"''. At the same time, in 
the regions of reduced number density between the spark- 
associated columns of plasma, where k ~ 10-^100, we obtain 
tp ~ 0.1. For the particle to emit the synchrotron radiation 
the condition 



(7) 



should be satisfied. Apparently, the latter condition is not 
met inside the dense plasma {ipf^ ~ 1 for 7^, ~ 100), whereas 
in the space of reduced number density it is satisfied well. 
In other words, only the particles, being in the space be- 
tween the spark-associated plasma columns, obtain (due to 
cyclotron absorption of radio waves) pitch-angles big enough 
to start gyrating relativistically, hence, emitting synchrotron 
radiation. 



3 THE SYNCHROTRON MODEL 

The total radiated power of a single electron synchrotron 
radiation is defined as (Ginzburg 1979; Lang 1980) 

Lo ~ 1.6 X 10"^^_B^7p sin^ V [erg s"^] , (8) 

and the critical frequency of the synchrotron radiation is 

i^c ~ 4.2 X 10^57^ sin 7/; [Hz] . (9) 
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The spectral density of a single electron radiation power near 
i^c is 

Lo 



lo (j^c) 



3.8 X 10~^^BsinV' [ergs"^Hz"^] 



(10) 



Let us notice that, if a source as a whole is moving towards 
an observer, the power given by equation (|) should be di- 
vided by sin^ ip (see equation (5.12) in Ginzburg (1979) and 
the discussion on page 75). However, this is not the case in 
our current application, as the emitting volume represented 
in Fig. ^ is quasi-stationary with respect to an observer on 
Earth. Indeed, the particles are permanently flowing in and 
out but the volume as the whole does not approach an ob- 
server with the relativistic velocity. That is why the multi- 
plier 1 / sin^ ip should be omitted in this consideration. 

The average magnetic field in the region where syn- 
chrotron emission is generated (Fig. |l]) can be estimated 
from the equation 



B = 3.2 X 10' 



(PP) 



a \ r) 



9.4 P- 



TZ^ sin a 



(11) 



as B ^ 3.3 X 10^ G. This provides the critical frequency 
i/c ~ 1.4 X lO'^^ Hz (equation which falls into infrared 
band (the frequency corresponding to the maximum in the 
single electron synchrotron spectrum yields i/m ~ 0.29 i/c ~ 
4 X 10^^ Hz). From equation (|l|) we find that Jo (uc) « 
1.3 X 10~^^ erg s~^ Hz~^. It is natural to assume that the 
spectral density of total power radiated by A*' electrons is 



/ (i/) = NIo (u) . 



(12) 



Let us then estimate A'^, i.e., the total number of emitting 
particles. From Fig. Q assuming the conservation of the 
particle fiux along the dipolar fields lines, we find that in 
the cut-cone with the height I = (7^2 — TZi) P^c cospi — 
1.3 Plc cos pi « 1.3x10^ cm there are TV KiZ ^ 4xl0^"«: 
particles (where Fpi is defined from equation p|) at any given 
moment of time. Then, from equation (^2|), for k — 10, we 
obtain 



I{uc) = 6 X 10^^ erg s"^ Hz' 



(13) 



which is a theoretical value of spectral density near the crit- 
ical infrared frequency Vc- The spectral density of power in 
the P-band can be deduced using the power-law spectrum 
of Geminga's optical emission (Martin et al. 1998), as 



3 X 10" erg s"^ Hz"\ (14) 



where = 6.9 x lO" Hz. 

Let us compare the theoretical value of I {Vg ) (equa- 
tion [l^ ) with that provided by the observations. The latter 
can be derived from the spectral density of fiux f {u^) = 
2 X 10"^° erg s"^ Hz"^ cm"^ (see Fig. 3 in Martin et al. 
1998), which gives the spectral density of total power 



"^fi^s) erg s ^ Hz \ 



(15) 



where ? is a beaming factor and d ~ 5 x 10'^'^ cm. According 
to Fig. ^ using the features of the dipolar geometry, we 
find that (f> ~ 1.5 pi ~ 0.8 rad, which leads to the beaming 
factor <; ^ (fp ^ 0.6 sr (instead of 47r sr for the isotropic 
radiation). Therefore, equation (|l|) yields I°'°''{Vs) « 3 x 
10^^ erg s~^ Hz~^, which is consistent with our theoretical 
prediction (equation ^) . 




Radio emission 
generation region 



Figure 2. Geometry of Geminga's optical radiation observed on 
Earth, f) is an impact angle of the closest approach of the ob- 
server's line of sight to the magnetic axis y,. Low-frequency radio 
emission which is not damped in the magnetsopheric plasma is 
also marked. See explanations in the text. 



Obviously, the fraction of the synchrotron optical radi- 
ation of Geminga observed on Earth is emitted in the same 
direction as the low- frequency radio waves (Paper I). Us- 
ing the feature of dipolar field lines that all of them in- 
tersect any given radial line with the same angle, we find 
a geometrical place of all the points in the magnetosphere 
where the tangents to the field lines point to the same di- 
rection^ as a wavevector fcr of the 100 MHz radio waves. 
As it is seen in Fig. ^, these points are placed at the radial 
straight lines inclined to /i at an angle « 13° from the 
side of the nearest approach to /i and ■& ^ 20° from the 
furthest one. The 'optically bright' region is restricted in 
the longitudinal direction to the altitudes along these lines, 
where the cyclotron damping of 0.1 -r- 1 GHz radio waves 
occurs. Calculations show (Khechinashvili et al. 1999) that 
the longitudinal size of this region ~ P^^ . The dimension 
of this region in the transverse direction can be estimated 
as l± ^ Rc/^p ^ 10~^Plc, where Pc is a magnetic field 
line curvature at corresponding altitude. Such an estima- 
tion of l± assumes that this is a characteristic distance at 
which a particle moving along the curved magnetic field line 
keeps emitting towards the observer (i.e., an observer stays 
in the emission diagram of a singel particle, within the an- 
gle range ~ l/7p). The corresponding characteristic time 



is Tx ~ Zj_/c ~ 4 x 10~ s, during which a particle radiates 
Eo = L(iT±_ energy. Obviously, the latter should be much less 
than the kinetic energy of an electron = 'y^mc^. Let us 
notice that in such a consideration the multiplier 1/sin^?/) 
(see the discussion below equation [lo| ) should be taken into 
account in equation as the source of the synchrotron 
radiation (a particle) is moving towards the observer in this 
case. Using equation (^) with the latter correction we find 



that X = Eo/Ek 



10" 



^ 1, so that only a tiny frac- 



t /3 = 20° at the nearest approach of the line of sight to the 
magnetic axis fj,, and /3 -|- 2o = 30° — at the furthest one. 
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tion of the particle's kinetic energy is radiated away as a 
synchrotron emission. Assuming that this is valid for all the 
radiating particles, we can estimate the total radiated power 
of the synchrotron radiation as L = xLpi ~ 8 x 10^^ erg s~^ 
(where Lpi is found from equation Bl). Then, using equa- 
tion (Q, we obtain I [u^) ^ 6 x lO^^rg s"^ Hz"^ We see 
that this value is consistent with our previous theoretical 
estimation (equation and hence with the observation- 
ally derived value (equation ^|). The consideration based on 
the energy conservation provides an independent test of our 
model. 

In the calculations above we assume that both pitch- 
angles tp and Lorentz factors 7^ of emitting particles remain 
constant while the particles cross the 'optically bright' re- 
gion. Obviously, this assumption is valid for Lorentz factors 
because x ^ 1- As for pitch-angles, it requires that a single 
particle goes through at least few acts of radio-absorption 
during the time Tm = l\\/c ~ 0.04 s. This is equivalent to 
fulfillment of the following condition 

Td/rm < 1. (16) 

Here = l/F is a characteristic time-scale of the cyclotron 
damping, and F is the damping decrement. We found (Pa- 
per I; Khechinashvih et al. 1999) that FR/c ~ 100 at the dis- 
tances where the bulk of radio band (0.1-^1 GHz) is typically 
damped. Therefore, F ~ 2 x 10^ s"^ and ra/r^ ~ 10"^ < 1. 
Thus, the parameters 7^ and 1/' used in our model can be 
considered quasi-stationary. 

We believe that the apparent optical spectral index 
a ~ 0.8 (Martin et al. 1998) indicates a power-law distri- 
bution of emitting particles over energy n (7) ~ 7 ', where 
g = 14- 2a f» 2.6 (Ginzburg 1979; Lang 1980). Such a dis- 
tribution function can be expected in the region where the 
synchrotron optical radiation of the Geminga pulsar is be- 
lieved to originate. 
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4 CONCLUSIONS 

The model for the optical radiation of the Geminga pul- 
sar, presented in this paper, is a direct consequence of our 
model for its erratic radio emission (Paper 1). Let us sum- 
marise our main results: (i) Relativistic charged particles of 
the pulsar magnetosphere obtain non-zero pitch-angles due 
to cyclotron damping of radio waves; (ii) The particles re- 
emit the absorbed radio energy in the form of synchrotron 
radiation, with a maximum power in the infrared band; (iii) 
The observed optical emission of the Geminga pulsar is just a 
short-wavelength continuation of its power-law synchrotron 
spectrum; (iv) The optical power estimated from our model, 
using two different approaches, agrees with the observed 
value; (v) We predict a significant pulsed infrared emission 
which should be observed from the Geminga pulsar; (vi) As 
a by-result, we support a model of the non-uniform plasma 
outflow in pulsar magnetospheres, originating in the form of 
spark discharges above the polar cap. 
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